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Abstract 

The effects of the symmetry energy softening of the relativistic mean field (RMF) models on the 
properties of matter with neutrino trapping are investigated. It is found that the effects are less 
significant than those in the case without neutrino trapping. The weak dependence of the equation 
of state on the symmetry energy is shown as the main reason of this finding. Using different 
RMF models the dynamical instabilities of uniform matters, with and without neutrino trapping, 
have been also studied. The interplay between the dominant contribution of the variation of matter 
composition and the role of effective masses of mesons and nucleons leads to higher critical densities 
for matter with neutrino trapping. Furthermore, the predicted critical density is insensitive to 
the number of trapped neutrinos as well as to the RMF model used in the investigation. It is 
also found that additional nonlinear terms in the Horowitz-Piekarewicz and Furnstahl-Serot-Tang 
models prevent another kind of instability, which occurs at relatively high densities. The reason is 
that the effective a meson mass in their models increases as a function of the matter density. 

PACS numbers: 13.15.+g, 25.30.Pt, 97.60.Jd 
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I. INTRODUCTION 



Recently, the dynamical stability of the uniform ground state of multi components sys- 
tems Q.e^ electrons, protons, and neutrons) at low densities has received considerable atten- 
tions taSQfly. This interest is motivated by the fact that the .eut.on St. is expected 
to have a solid inner crust of nonuniform neutron-rich matter above its liquid mantle ^ and 
the mass of its crust depends sensitively on the density of its inner edge and on its equation 
of state (EOS) [2J. Meanwhile, the critical density (pc), a density at which the uniform 
liquid becomes unstable to a small density fluctuation, can be used as a good approximation 
of the edge density of the crust 0. Using the Skyrme SLy effective interactions. Ref. Q 
found the inner edge of the crust density to be pedge = 0.08 fm ^. Reference j3| generalized 
the dynamical stability analysis of Ref. in order to accommodate the various nonlinear 
terms in the relativistic mean field (RMF) model of Horowitz-Piekarewicz 0| . They found 
a strong correlation between pc in neutron star and the symmetry energy (agym) which leads 
to a linear relation between pc and the skin thickness of finite nuclei. 

These results suggested that a measurement of the neutron radius in ^osp]^ ^^]^ provide 
useful information on the Pc Ql- Recently, a new version of the RMF model based on 
effective field theory (ERMF) has been propo.ed by Fnrn.ah,. Serot. arrd Tan, g The 
predictive power of the model in a wide range of densities is quite impressive (see the review 
articles 0, Q| for the details). On the other hand, an adjustment of the isovector- vector 
channel of the ERMF model in order to achieve a softer density dependence of Osym at high 
densities has been done in Ref. 

Similar to the previous calculation done in P], but using a different RMF model, it is 
found that by adjusting that channel a lower proton fraction (1^) in high densit y n eutron 
star matter can be achieved without significantly changing the bulk properties It is 

well known that Yp is related to the threshold of the direct URCA cooling process and 
the trend of the density distribution of Yp is unique for each model which is sensitive to 
the different forms of nonlinear terms used. Therefore, an investigation of the EOS and 
the instability at low densities by using different RMF models could be quite interesting. 
Indeed, by comparing our results with the previous ones [sl, 0] , we can systematically study 
the influence of the form and the strength of isovector-vector nonlinear terms on the critical 
density. The presence of electrons and the influence of the electromagnetic interaction on 
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the unstable modes of asymmetric matter as well as the comparison between dynamical 
instability region with the thermodynamical one have been investigated in Refs. 3 0] by 
using the standard RMF model within the Vlasov formalism. They observed the important 
role of the Coulomb field in the large structure formation and also the role of the electron 
d.„a,„,csi„.eston.gthe.a.ge-wavele.gthi„.tab„,tiesa. 

In the actual situation, muons may also be present in the stellar matters. Their existence 
yields an additional electromagnetic contribution besides the electrons and protons ones. 
The significance of their contribution depends on the matter composition. Including muons 
in the dynamical instability analysis of uniform matter increases the electromagnetic effect on 
the Pc in medium. Moreover, in supernovae or protoneutron stars neutrinos can be trapped 
inside, if their mean free paths are smaller than the star radius. In general, the presence 
of trapped neutrinos in matter affects the stiffness of EOS and drastically changes the 



composition of neutral matter ll^ 
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15, m. It was a 



so found that neutrino trapping 
to appear at higher densities 
, the EOS of the strangeness 



causes the strange baryon [3, Il3| and kaon condensate [3, 
as compared to the case without neutrino trapping. In Ref. 
rich protoneutron star and the EOS of neutron star matter with temperature different from 
zero have been already studied by means of the ERMF model. Meanwhile, like neutron stars, 
protoneutron stars have dense cores and outer layers Q| and in supernovae inhomogeneity 
can appear below the saturation density of nuclear matter (po) jl?! - Therefore, an 

investigation of the boundary between the two phases, which leads to the determination of 
Pc becomes crucial for a realistic and complete description of stellar matters. 

In thepresent work, we shall extend the analysis of dynamical stability in uniform matter 
of Ref. p in two ways, i.e., first, we consider the muon contribution and, second, we shall 
consider all possible mixings of vector and isovector contributions caused by the nonlinear 
terms in the ERMF model. The parameter set NL3 of the standard RMF Il9|l model as well 

n 

as the parameter sets Z271 and Z271* of the Horowitz-Piekarewicz one |J| will be presented 
for comparison. The results are used to study the effect of the neutrino presence on the pc 
values. Here, we only use nucleons in the baryonic sector along with a, cu, and p mesons in 
the mesonic sector because pc usually appears at a density lower than po and the strange 
particles appear at a density higher than 2po. Furthermore, in the neutrino trapping case 
they appear at higher density than in the case without neutrino trapping. Therefore, we 
can assume the strange particles to have only minor effects on the pc and, as a consequence. 
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their contributions may be neglected in this work. On the other hand, although in the real 
situation the temperature of protoneutron stars is not equal to zero and the supernovae inner 
core can have a temperature around T ~ (10-50) MeV, the zero temperature approximation 
can be assumed here since the temperature effect on the EOS of supernovae matter and 



on the maximum mass of protoneutron star 



14 Il5| is smaller than that without neutrino 



trapping. In this approximation, the following constraints can be used to calculate the 
fraction of every constituent in matter: 

• balance equation for the chemical potentials 

/^n + l^u, = /^p + l^e, (1) 

• conservation of the charge neutrality 

Pe + Pt, = Pp, (2) 

where the total density of baryon is given by 

Pb = Pn + Pp, (3) 

while the fixed electronic-leptonic fraction Yi^ is defined as 

Y,^ = ^^^^Y, + Y^^, (4) 
Pb 

where Y^^ and Yg are neutrino electron and electron fractions, respectively. In addition, we 
will revisit the EOS of matter with neutrino trapping because we want to study the effect of 
different adjustments in the isovector- vector sector [llj on the 
trapping. They were not fully explored in our previous works UJ, |2i| . 

This paper is organized as follows. In Sec. II, a brief review of the model is given. 
Calculation of the dielectric function is presented in Sec. Ill, while numerical results along 
with the corresponding discussions are given in Sec. IV. Finally, we give the conclusion in 
Sec. V. 
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II. RELATIVISTIC MEAN FIELD MODELS 



Our starting point to describe the non-strange dense stellar matter is an effective La- 
grangian density for interacting nucleons, a, u, p mesons as well as noninteracting leptons. 

4 



The Lagrangian density for standard RMF models can be found in Refs 



.0 
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m. We 



note that in order to modify the density dependence of agym, Horowitz-Piekarewicz have 



recently added isovector- vector nonlinear terms (£hp) in a standard RMF model 
case of ERMF model, the corresponding Lagrangian density can be found in Refs. 



. In the 
11, d. 



For the convenience of the reader, we write the Lagrangian density in a compact form as 



C — Cn + + -^HP + ^L, 

where the nucleon part, up to order z/ = 3, reads 



(5) 



with 



P 



n 



(6) 



(7) 



(8) 



^ = exp(27r(x)//^), 7i{x) = -f ■ n{x), 



(9) 



(10) 



(11) 



(12) 



(13) 



where, p and n are the proton and neutron fields, M is the nucleon mass, while a, vr, V^, 
and are the cr, tt, lj and p meson fields, respectively. The meson contribution, up to order 
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■•M 



\flTi{d,Ud^U^) + \flTi{UU^ - 2) + \d,ad^a 
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+ ( 1 + j mjTr(6^6'^) - m'y [ 1 + 



3! M 4! M2 



where 



f/ = 



(14) 



(15) 



In the mean field approximation, vr meson does not contribute. In the Horowitz-Piekarewicz 
model 41 the isovector-vector nonlinear term reads 



£hp = ^^vqIqI ■ b, V^V,. 



(16) 



We note that changmg the and A. parameters m the Horow.tz-P.ekarew.z model 
or gp and rip parameters in the ERMF model [111], affects the density dependent part of the 
. A more detailed procedure to adjust the density dependence of the nuclear matter 



•^sym 



symmetry energy in RMF models can be found in Refs. j4, 
For the leptons, the free Lagrangian density 



11 



24|. 



l=e- , ^- , u^, 

is used. In the present work, the following parameter sets are chosen: 

• Standard RMF model: NL3 Q, 

• Standard RMF model plus isovector-vector nonlinear term: Z271 and Z271* jj], 

• ERMF models: G2 M and G2* \ 



(17) 



The coupling constants for all parameter sets are shown in Table lU 
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TABLE I: Numerical values of the coupling constants used in the parameter sets. 



Parameter 


G2 


NL3 


G2* 


Z271 


Z271* 




0.554 


0.541 


0.554 


0.495 


0.495 


9a/(4^) 


0.835 


0.813 


0.835 


0.560 


0.560 


5^/(4^) 


1.016 


1.024 


1.016 


0.670 


0.670 


5p/(4vr) 


0.755 


0.712 


0.938 


0.792 


0.916 




3.247 


1.465 


3.247 


1.325 


1.325 




0.632 


-5.668 


0.632 


31.522 


31.522 


Co 


2.642 





2.642 


4.241 


4.241 




0.650 





0.650 








V2 


0.110 





0.110 










0.390 





4.490 








Ay 














0.03 



III. DYNAMICAL INSTABILITY OF UNIFORM DENSE STELLAR MATTERS 



In order to accommodate muons and various nonlinear terms in the ERMF model, we 
have to extend the stability analyses of the uniform ground state of Refs. 0, Q, 22 
end, the longitudinal polarization matrix in these analyses is modified as follows 



To that 



/n 



00 







n TT^ 

u ■'^-'■00 



TT W TT" 



(18) 



\ 

_ 



n^o 

V n;: w^j 

where the scalar polarization lis = 11^ + 11^ . In Eq- (|18|)- the polarization due to the mixing 
between scalar and vector terms through sigma meson is indicated by 11^. In the limit of 
0, the individual polarization components are 



ns(g,o) 



+ 



27r2 ^ 



kpEp 
In 



3M* 2 + — 



In 



kp + Ep 



2kF — q 


2E\ 
m 


2kF + q 


q 



M* 

qEp - 2kpE 



qEp + 2kpE 



(19) 
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for the scalar polarization, 



n^(g,o) 



M* 
2^ 



In 



2kf 



1 



2kF + q 



(20) 



for the mixed scalar-vector polarization, and 



noo(g,o) 



1 



2, p q\ kp + Ep 

-kpEp In 

3 6 M* 



El 

3q 



M* 2 + 4 - 



3g2 
4 



In 



2kf 



q 



2kF + q 



E 



q' 



+ _ ( M* ^ - ^ 1 In 



qEp - 2kFE 



(21) 



qEp + 2kFE 

for the longitudinal polarization, with Fermi momentum kp, nucleon effective mass M* = 
M - g^a, Fermi energy Ep = {kj, + M* 2)i/2 ^nd E = {q^/A + M* 2)V2_ por electron (e) or 
muon (/i), M* is equal to electron or muon mass, respectively. 
The longitudinal meson propagator now reads 

[ dg dg -dg \ 

dg dg —dg 

Dr = -4 dl^ d-^p , (22) 
~-dg —dg d~^„„ d. 
\ d: 



vp 



svp 



d^p (^44 y 



where = -{dsv + dsp), dj^^ = -{dsv - d^p), d;;^ = d^ - dp, d^s = dg + d^ + dp + 2d^p and 
dii = d^ + dp — 2dyp. In this form, mixing propagators between isoscalar-scalar and isoscalar- 
vector {dgy), isoscalar- vector and isovector- vector (d^p), isoscalar-scalar and isovector- vector 
(dsp) are present due to the nonlinear terms in the model, in addition to the standard photon, 
omega, sigma and rho propagators {dg, d^, dg and dp). These propagators are determined 
from the quadratic fluctuations around the static solutions which are generated by the second 
derivatives of energy density {d'^e/ d(j)id(f)j), where (pi and (pj are the involved meson fields. 
The energy density e derived from Eq. (0) reads 



(2vr) 



^ j d^kiEi{ki) + gu,Vo{pp + pn) + ^Qpboipp 



Pn) 



t=p,n,e,fi,i/e,i^ij. 
1 



2 4 



(23) 
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From Eq. ()23p we can obtain the explicit form of all contributions to the longitudinal 
propagator. Note that the energy density of the ERMF model can be obtained from Eq. 
by using the following explicit expressions of the coupling constants: 



2M ' ^ 6M2 ' 2M 





A. = A. = 0. (24) 

On the other hand, if we set the coupling constants in Eq. ()23|) to 

f^ = d2 = ds = 0, A, = 2A.,glgl, A, = 2A,glgl (25) 

we will obtain the energy density of the Horowitz-Piekarewicz model . The mesons effective 
masses and the mixing polarizations calculated from the energy density [Eq. (j23p ] are given 
by 

ml ' = = ml + 2d2a + d^a^ + ScsFq' + A.&o. 



(26) 



and 



= -a^ = 2^* + ^^"^'»• 



By substituting Eqs. ()26jl and (j27|l in the a, u, and p propagators, 

^ ^ t/g(g^ + m* ^)(g^ + m; ^) 



(27) 



(28) 



- (g2 + ^*2)(^2 + ^*2)(^2 + ^*2) + (n0,)2(^2 + ^*2)_(n00)2(g2+,„*2)' ^ ^ 
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, ^ ll^g%q' + m*J){q^ + m*J) 



and in the mixing propagators, 



i/(g,go = 0) 



m 



* 2> 



with 



i7(g,go = 0) ' ^^^^ 



+ iKfici' + <')-iT^Znci' + K% (34) 



and photon's propagator 



c/g = ^, (35) 



and using the exphcit form of each component, the longitudinal meson propagator can be 
obtained. 

The explicit derivation of each propagator is given in Appendix ^ Note that by setting 
all coupling constants in Eqs. (j^Hj) - ()33|) which are not required by the Horowitz-Piekarewicz 
model 0] to zero, Eqs. (8), (13) - (14) of Ref. can be obtained. 

The uniform ground state system becomes unstable to small-amplitude density fluctua- 
tions with momentum transfer q when the following condition is satisfied 

det[l-Di(g)ni(g,go = 0)] <0. (36) 

The explicit form of Eq. ()3(ij) is given in Appendix El In the case that the density is smaller 
than pq, the critical density pc is the largest density for which Eq. (jHBj) has a solution. In 
the case that the density is larger than po, if any, p^ is the smallest density. 
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FIG. 1: Binding energy per nucleon (upper left panel) and the neutrino electron fraction (upper 
right panel) as a function of the baryon density, as well as the equation of states (EOS) as functions 
of baryon density (lower left panel) and energy density (lower right panel) according to RMF 
models. The curves are obtained by using Yi^ = 0.3. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

Before discussing the dynamical instability of non-strange dense stellar matter, we will 
discuss the effects of different treatments in the isovector-vector sector on the EOS, 
binding energies (E/A) and the neutrino electron fraction (Y^^) of matter with neutrino 
trapping. The results are shown in Fig. ^ In contrast to the case of matter without 
neutrino trapping significant differences in the trends of E/A and EOS start to appear 
at Pb/po ~ 2.0. Except at sufficiently high densities, the difference between G2 and G2*, 
or Z271 and Z271*, does not significantly show up in these trends. These results can be 
understood from the asymmetry expansion of the binding energy of the corresponding matter 
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in the vicinity of the symmetric nuclear matter (SNM). The latter is given by 



25 



E/A{pB, a) = {E/A)s^^{pB, a) + a^asymips) + 0{a^) +AEl/A{pb), (37) 

AEi/A(pB,a) AE2/A{pB,a) 

where a = — Yp and AEl/A is the electrons and muons contribution to the binding 
energy. We have found that for all parameter sets and matters used the value of AEl/A is 
smaller than the other terms. The AE2/A term is given by 3| 

AE2/A{pB, a) = a^QipB) + ... (38) 

The origin and connection of the quartic term [Q{pb)] to direct URCA are discussed in 
Ref. |26|. For the case of pure neutron matter {AEi/A = asym) or other fixed a cases it is 
known that AE2 <^ AEi. This means that the convergence of the expansion is considerably 
fast for these asymmetric matters but for matter with and without neutrino trapping the 
situation is quite different (see Fig. We can see that by imposing the neutrality and jJ 
stability conditions on the matter with and without neutrino trapping, the corresponding 
asymmetry (a) becomes density dependent and, evidently, their AE2/A are substantially 
larger compared to that of the PNM, for example. This indicates that in these cases, the 
convergence is significantly slow or even can not be reached at all. Additional constraint in 
the form of fixed electronic lepton fraction (Yj^) for the case of matter with neutrino trapping 
causes the decrease of and the convergence of binding energy expansion are slower than 
those in neutrinoless matter. 

In Fig. El we show the AE2/A and AEi/A as a function of the ratio between baryon and 
nuclear saturation densities for the matter with and without neutrino trapping where G2, 
G2*, Z271, Z271* and NL3 parameter sets are used. The difference between the two types 
of matters appears mainly at densities lower than 2po, where AE2/A is smaller than AEi/A 
for neutrinoless matter. The opposite situation happens for the neutrino trapping case. In 
the case of neutrinoless matter and parameter sets with a stiff agym, AE2/A is larger than 
AEi/A for the density larger than 2po, but for the case of soft asym this condition is reached 
only after the density becomes relatively high. Figure E] shows the characteristic of for 
both matters as the baryon density increases. Here, the dependence of on parameter 
sets is obvious for matter without neutrino trapping, while the opposite situation happens 
for matter with neutrino trapping. The symmetry energy is shown in Fig. [SJ where we can 
clearly see its dependence on the parameter sets. The dependency of AE2/A, AEi/A and 
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FIG. 2: /S.E2/A, A£'i/A and as a function of the ratio between baryon and nuclear saturation 
densities for matters with and without neutrino trapping (NT). The G2* parameter set and Yi^ = 
0.3 are used in obtaining these results. 
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FIG. 3: ^E^l A and AE'i/A as a function of the ratio between baryon and nuclear saturation 
densities for matters with and without neutrino trapping (NT), obtained by using G2, G2*, Z271, 
Z271* and NL3 parameter sets. AU curves are obtained with Yj^ = 0.3. 

o? on asym can be investigated by comparing Figs. El and El with Fig. O It is obvious that 
in the case of matter with neutrino trapping, correlations between IS.E2I IS.E\I A and o? 
with ttsym at densities lower than 2po can not be observed, but for the neutrinoless case 
the opposite situation appears. For densities larger than 2po, correlations between Ai?2M, 
Ai^i/A, as well as o? and a^^m for matter with neutrino trapping is weaker than those 
without neutrino trapping. 

Therefore, it is clear that the results displayed in Fig. ^ are caused by one reason: matter 
without neutrino trapping has a weak correlation with agym and especially at densities lower 
than 2/9o, the effect is more pronounced. This effect is mainly due to the behavior of the 
IS.E2I A contribution. For densities lower than 2po, the asymmetry between protons and 
neutrons, which is smaller compared to the case of neutrinoless matter, also induces a 
correlation between [\E\I A and agym, whereas the role of l\E\l A is suppressed. 
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FIG. 4: as a function of the ratio between baryon and nuclear saturation densities for matters 
with and without neutrino trapping (NT). The results are obtained by using G2, G2*, Z271, Z271*, 
NL3 parameter sets and Yi^ = 0.3. 




FIG. 5: Symmetry energy as a function of the ratio between baryon and nuclear saturation densities 
for G2, G2*, Z271, Z271* and NL3 parameter sets. 

In the upper right panel of Fig. Q Y^^ for different parameter sets are shown. It is 
interesting to note that NL3, Z271 and G2 parameter sets have a similar Yi,^ trend but Y,y^ 
of G2* and Z271* behaves differently. By comparing this figure with Fig. Elwe can conclude 
that the number of trapped neutrino correlates with the agym behavior of the models. This 
is due to the neutrino fraction which self consistently depends on the proton fraction (Yp) 
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FIG. 6: Same as Fig. ^ but here the Yi^ is varied and the G2* parameter set is used. 

through the neutrality and beta equilibrium conditions, while Yp has a correlation with flsym, 
although the effect is quite small for this kind of matter. 

To study the role of the neutrino number in the E/A and EOS using G2* parameter set, 
we show in the upper left panel of Fig. El the variation of E/A for some values of Yj^ and 
in the lower panels the variation of their EOS. It is found that the larger the number of 
neutrinos in matter, the smaller the value of E/A. At low density, the EOS of matter with 
neutrino trapping is stiffer than that without neutrino trapping. According to Ref. the 
reason is that the neutrino in matter shifts the threshold of muon production toward higher 
densities. However, for matter with neutrino trapping, the larger the number of neutrinos, 
the softer its EOS. At high densities the situation is opposite. Also we can see this from 
another point of view, i.e., by comparing the slope of dominant contributions of each case 
(bottom and center panels of Fig. [7j). At densities less than po, the slope of AEi/A for 
neutrinoless matter is smaller than that of AE2/A for matter with neutrino trapping. In 
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FIG. 7: /S.E2/A, IS.E1/A and a as a function of the ratio between baryon density and nuclear 
saturation density for matters with neutrino trapping (NT) obtained by using the G2* parameter 
set but with varied Yi^ . 
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FIG. 8: Critical densities for the Z271, Z271*, G2 and G2* parameter sets as a function of the 
neutrino fraction in matter. 



the latter, also in the same density range, even though it does not clearly visible, the slope 
becomes larger as Yi^ becomes larger. The situation is reversed once the densities become 
higher than pQ. 

The critical densities as a function of the neutrino fraction in matter are shown in Fig. |S1 
for the Z271 (upper left panel), Z271* (upper right panel), G2 (lower left panel), G2* (lower 
right panel) and in the left panel of Fig. El for the NL3 parameter set. Without neutrino 
trapping, it is obtained that p)?^ = o.052 fm"^ p^^* ^ 0.O6O fm~^ pf = 0.066 fm"^ 
^Z27i* ^ Q_Qg2 pNL3 ^ 049 f]^-3_ rjj^jg indicates that the critical density in the 

Horowitz-Piekarewicz model is larger than that of the ERMF and the NL3 parameter sets, 
while the standard RMF model represented by the NL3 parameter set yields the smallest 
value. The effects of the isovector-vector adjustment on pc can be seen by comparing the 
panels for Z271 with Z271*, or G2 with G2*. In the Horowitz-Piekarewicz model, the 
variation of the pc values is larger than in the ERMF one. 
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FIG. 9: Same as in Fig. |S1 but for the NL3 parameter set. 




FIG. 10: Effective masses of micleon, a, oo and p mesons as a function of the baryon density. Note 
that several Unes coincide. 
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On the other hand, with neutrino trapping we obtain, p^^ = (0.084 — 0.086) fm ^, 
pG2* ^ (0.083-0.086) fm-^ pf^^ = (0.085 - 0.086) fm-^ pf^i* = (0.087 - 0.088) fm"^ 
and p^^^ = (0.081 — 0.084) fm~^. This indicates that for this case the value of pc does not 
sensitively depend on the model and on the variation of the number of trapped neutrinos in 
matter. In general, their values are larger than those obtained in the case without neutrino 
trapping. 

Equation (|36p implies that pc is determined by matter composition as well as by the 
effective masses of mesons and nucleons used. We know from previous results that protons 
and neutrons are the dominant constituents in matter composition and its behavior can be 
observed from the corresponding asymmetry («) which has been discussed and shown in 
Figs, m and IHl On the other hand, the effective masses depend also on the used parameter 
sets (see Fig. ITO|l . The above findings and the facts that for each parameter set p^'**^ < 

pwithout NT^ 

well as the indication that in matter with neutrino trapping and constant Yi^ 
the critical density does not vary significantly for all parameter sets used, can be understood 
as the evidence of the dominant role played by the proton-neutron asymmetry in matter in 
determining the value of pc- Furthermore, the reason that the strong dependence of pc on 
the parameter set used in the case of matter without neutrino trapping is that the protons- 
neutrons asymmetry of this kind of matter has a strong correlation with the ctsym- 

From Figs. El Eland the top panel of Fig. [71 we can clearly see that the dependence of Pc 
on the number of trapped neutrinos in matter does not significantly change with respect to 
the variation of Yi^ and the value of a^. This result emphasizes our previous finding that 
there is a strong relation between p^ and the agym- 

To show the crucial role of protons in shifting the values of pc to relatively low densities, 
let us study the contribution of each constituent. First, the electrons, protons, and muons 
contributions are switched off. Then, by fixing q at 30 MeV we search for pc for all parameter 
sets, both for trapped and untrapped neutrino cases. For all cases, we found more or less 
similar results, i.e., Pc/po ~ 1-7. This value comes merely from the neutron contribution. 
A small difference in the neutron fraction {Yn) for both cases has a negligible effect on the 
Pc prediction. Secondly, now in addition to the first condition, the proton contribution is 
switched on. Now we obtain quite different values of pc, for example, in the case of G2* 
with Yi^ = 0.3 we obtain Pc/po = 0.580, and for the case without neutrino trapping we 
obtain pc/ po = 0.425. Thus, it is clear that the protons shift the pc values to the region with 
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relatively low densities. Next, the electron contribution is switched on. Then we observed 
that Pc/po is shifted further, i.e., Pc/po = 0.430 for Yj^ = 0.3 and Pc/po = 0.350 for matter 
without neutrino trapping. If the muons contribution is switched on, the result does not 
change. It happens not only in matter with neutrino trapping, but also in matter without 
neutrino trapping, respectively. However, these results show that the proton and electron 
contributions decrease the value of pc in different ways for both cases and, therefore, reveal 
the crucial role of matter composition in determining the value of pc- 

In contrast to the Horowitz-Piekarewicz and ERMF models, the NL3 parameter set yields 
an additional instability at relatively high densities. This is shown in the right panel of Fig. IHl 
for the g = 2 MeV case, i.e., pd Po ~ 6.5. Even if we use q close to zero, this instability does 
not disappear and the pd Po value stays similar as in the q = 2 MeV case. Therefore, this 
instability is not due to the small-amplitude density fluctuations. In fact, this instability 
appears because the effective a mass of the NL3 parameter set is zero at that density (as 
shown in the lower right panel of Fig. lTUI) . As a consequence, the corresponding a propagator 
changes the sign at this point. This fact clearly shows that additional nonlinear terms to the 
a nonlinearities of the standard RMF models like in the Horowitz-Piekarewicz and ERMF 
ones can prevent such kind of instability to appear. 

V. CONCLUSION 

The effects of the different treatments in the isovector-vector sector of RMF models on 
the properties of matter with neutrino trapping have been studied. The effects are less 
significant compared to those without neutrino trapping. Different dependences of the EOS 
and B/A of both matters on agym are the reason behind this. 

The effects of the variation of the neutrino fraction in matter on the EOS and B/A have 
been also discussed. 

The longitudinal dielectric function of the ERMF model for matter consisting of neutrons, 
protons, electrons, muons, and neutrinos has been derived. The result is used to study the 
dynamical instability of uniform matters at low densities. The behavior of the predicted 
Pc in matter with and without neutrino trapping has been investigated. It is found that 
different treatments in the isovector-vector sector of RMF models yield more substantial 
effects in matter without neutrino trapping rather than in matter with neutrino trapping. 
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Moreover, for matter with neutrino trapping, the value of pc does not significantly change 
with the variation of the models as well as with the variation of the neutrino fraction in 
matter. In this case, the value of pc is larger for matter with neutrino trapping. These 
are due to the interplay between the major role of matter composition and the role of the 
effective masses of mesons and nucleons. It is also found that the additional nonlinear terms 
of Horowitz-Piekarewicz and ERMF models prevent another instability at relatively high 
densities to appear. This can be traced back to the effective a mass which goes to zero when 
the density approaches 6.5 po • 
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APPENDIX A: MESON PROPAGATORS 

In this Appendix, simplified Dyson equations for meson propagators (only zero component 
of each propagator is considered) will be given. The covariant form for a and uj couplings is 
given in Ref. j28|. If we define the a, uj and p free propagators as 

r = ^ r °° = r °° = ~^ ( ^^\ 

the Dyson equation for the a propagator in the absence of coupling to uj and p fields is 
obtained by considering the sum of ring diagrams, i.e., 

G„ = G„- iG„Ii„„G„, (A2) 

from which we can obtain 

Similarly, the Dyson equations for the zero components of uj and p meson propagators are 

-i 
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and 



G^^ —i 



where {m*)^= {m^f +U^^, {mlf= (m^)^ +UZ, and (m;)2= (m^)^ +nOO. 

In the ERMF model, each meson is coupled to other mesons through the nonlinear terms. 
This fact further complicates the form of the full propagators. The Dyson equation for the 
a propagator including the possible mixing terms becomes 

Ga = — G^U^^^G^U^^^G^ — Go-n°p(5°°n°^Go-, (A6) 

which can be written as 

i + G.(nLGoonL + no^Goono^)- ^ ^ 

Equation ()A7|) can be simplified into 

i(^q^ — m*^ ^)(?^ — ^* ^) 
^ iQl - ml ^){ql - m; 2)(g2 _ 2) + (n^ J2(gf - ^) + (no^)2(g2 _ 2) ' (^8) 

Similarly, the zero component of uj and p meson propagators can be written as 

^00^ 

- (g2 - ^* 2)(^2 _ ^* 2)(^2 _ ^* 2) + (no J2(^2 _ 2) _ (n00)2(g2 _ ^* 2) ' ^ ^ 

and 

^00 ^ .A. 01 

" (g^ - m* 2)(g2 _ 2) (^2 _ 2) + (n0^)2(g2 - 2) _ (n00)2(g2 _ 2) ' I ^ 

We can also define a propagator GJJ,^ which contains the sum of all diagrams that transform 
UJ into 0", i.e., 

+ ■ ■ .. (All) 

This propagator may be summed up to produce 

_,-r'OOnO n 

/-yO _ '''-^1^ ^^uja^<T ('A12') 



which can be simplified to 
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Similarly, we can also obtain the sum of all diagrams that transform p into a as 

and that transform p into uo as 

where 



Equations (|1H1) - (|SS1) are special cases of Eqs. (IXsl- 1X101 rnairnHMBll . i.e. by taking the 
limit of go ^ and inserting the proper mesons coupling constants in the latter. 



APPENDIX B: EXPLICIT FORM OF THE LONGITUDINAL DIELECTRIC 
FUNCTION 



In this Appendix, the explicit form of the longitudinal dielectric function 
[1 — DL{q)IlL{q, qq = 0)] is provided. If we define the matrix ej^ as 

f Ai Bi Di El Fi \ 
A2 B2 D2 E2 F2 

A4 i?4 D/^ £^4 F4 

A5 -B5 -D5 F5 
V ^6 De Ee / 

then the contents of each component of the matrix in Eq. (jBljl are 



Ai 


— 1 _ // 


A2 


— — rl 


A, 


= 0, 


Bi 




B2 


= 1 — dgll^Q, 


B, 


= 0, 


Di 


— (] W 


D2 


— d W 




— 1 _L /7 n — d~^ w 


El 


— UgllQQ, 


E2 








Fi 


= 0, 


F2 


= 0, 


F4 


— // TT" — (7~ TT" 



(BI) 



d~ U"^ 

^svp my 



(B2) 
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and 



E,^i-d+^w^-dssni„ 

F5 = -dtvp^m - 



Ae = 0, 

^6 = 0, 

= -djvp^s - d~pW^ - duU"^, 
Ffi = 1 - dZ,Jll - duU^So- 



(B3) 



svp m 
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